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CURRENT  ENHANCEMENT  IN  A  CONDUCTING  CHANNEL 


In  a  series  of  recent  experiments , 1 it  was  observed  that  when  a 

relativistic  electron  beam  propagates  in  a  gas  (air,  neon,  or  nitrogen)  at 

about  100  torr,  the  net  current  is  higher  than,  and  at  times  almost  twice  as 

large  as  the  beam  current.  This  is  remarkable  as  the  plasma  current  generated 

by  the  beam  ionization  usually  flows  opposite  to  the  beam.  Concurrent  with 

this  current  enhancement  is  detected  a  large  amplitude  hose  motion  of  the 

electron  beam.  Several  attempts  have  been  made  to  explain  this  current 
14  R 

enhancement.  ’  Our  recent  numerical  simulations3  show  that  large  amplitude 
hose  motion  can  enhance  the  current  up  to  eighty  per  cent,  in  reasonable 
agreement  with  experimental  observations . 1 

In  this  paper,  we  show  that  current  enhancement  is  a  general  feature  when 
a  current  source  (electron  beam)  moves  sideways  inside  an  electrically 
conducting  channel.  Specifically,  the  present  study  indicates  that  the  total 
current  (the  beam  current  and  the  induced  current  in  the  plasma)  may  be  as 
much  as  three  times  the  .  beam  current  when  the  beam  reaches  the  conducting 
boundaries  of  the  channel.  These  conclusions  are  valid  as  long  as  the  channel 
conductivity  is  sufficiently  high.  They  are  a  direct  consequence  of  the 
Ampere-Faraday 's  Laws  and  appear  to  be  relatively  insensitive  to  the  geometry 
of  the  channel.  Following  the  presentation  of  our  analysis  of  a  highly 
simplified  model,  we  provide  an  interpretation  of  this  factor  of  three  of 
current  enhancement.  The  implications  will  be  addressed  but  the  details  will 
be  published  elsewhere. 


To  illustrate  the  ideas,  consider  first  a  one-dimensional  system.  In  our 
model,  an  infinitesimally  thin  current  sheet,  of  surface  density  Ig,  flows  in 
the  z-direction.  This  current  sheet  is  initially  located  midway  between  two 
perfectly  conducting  plates,  Pg  and  PR,  at  x  =  0  and  x  =  L,  respectively.  We 
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assume  that  the  medium  between  and  Pg  is  characterized  by  a  constant 
electrical  conductivity  o.  We  further  assume  that  the  current  sheet  has  been 
situated  at  x  =  xg  =  L/2  for  a  long  time  so  that  a  steady  state  is  already 

established  at  time  t  =  0,  i.e.,  any  plasma  current  has  decayed  to  zero  at  t  = 

0.  For  t  >  0,  the  current  sheet  is  given  a  uniform  x-ward  motion,  i.e.,  xg(t) 

=*  L/2  +  ut  where  the  speed  u  is  assumed  to  be  constant.  In  the  following,  we 

only  consider  the  response  of  the  magnetic  fields  to  this  transverse  motion  of 
Ig.  We  ignore  the  reaction  on  Ig,  due  to  the  resulting  modification  of  the 
field  profiles.  Presumably  this  sideward  motion  is  caused  by  some  beam 
instabilities,  or  simply  externally  driven. 

a 

In  terms  of  the  vector  potential  A  =  zA(x,t),  the  above  system  is 
governed  by  the  diffusion  equation 


c 

4it 


92A 

3x2 


-Iff*  V 


0<x<L,  0<t<  L/2u 


(1) 


where  c  is  tne  speed  of  light  and  5  is  the  Dirac  delta  function.  Gaussian 
units  are  used  and  displacement  currents  are  ignored.  The  boundary  condition 
to  Eq.  (1)  is  A(0 , t )  -  A(L,t)  =  0,  since  PL  and  PR  are  perfect  conductors.  The 
initial  condition  to  (1)  is  A(x,0)  =  (irLIg/c)  (1  -  j2x  -  L|/L)  for  0  £  x  £  L, 
reflecting  the  assumption  that  a  steady  magnetic  field  is  established  prior  to 
t  =  0. 

In  terms  of  the  vector  potential  A,  the  (inductive)  electric  field 

E  =  -(c  )  3A/3t  and  the  magnetic  field  B  =  V  x  A.  Thus,  the  term 

-(o/c)3A/3t  in  Eq.  (1)  is  the  current  density  induced  within  the  conductivity 

channel  and  the  current  enhancement  factor  F  is  given  by  F(t)  =1  - 
L 

(o/I^c)  J  dx  3A/3t,  which  is  the  total  current  within  the  channel  divided  by 
o 

the  beam  current. 
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The  diffusion  equation  (1),  together  with  its  associated  initial 
conditions  and  boundary  conditions,  may  be  solved  by  standard  techniques.  It 
can  be  shown  that,  in  the  limit  t  ♦  t  =  L/2u,  (i.e.,  shortly  before  the 
current  sheet  makes  contact  with  PR)  the  current  enhancement  factor  F  is  given 
by 


1  +  tanh  7  (-1 ) 

n=odd 


n-1 
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2 

n  tt 
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n(n£ 


v2) 


(2) 


where 


v  a  a  4uoL/c2  (3) 

2  2 

is  the  ratio  of  the  characterist ic  diffusion  time  t.  -  2oL  /c  and  the  transit 

d 

time  it  =  L/2u.  The  infinite  sum  in  Eq.  (2),  converges  rapidly  for  finite 
values  of  v. 

The  enhancement  factor  F,  as  given  by  Eq.  (2),  is  shown  in  Fig.  1  as  a 
function  of  the  normalized  velocity  v.  Several  features  are  noteworthy. 
Perhaps  the  most  interesting  one  is  that  F  approaches  an  asymptotic  limit  of 
three  as  v  becomes  large.  From  Fig.  1,  one  sees  that  F  >  2  for  v  >  1,  Thus, 
according  to  the  present  model,  the  total  current  within  the  channel  exceeds 
twice  the  beam  current  the  moment  the  current  sheet  strikes  the  plate  PR  if 
>  xt.  For  v  =  Tlj/,Tt  —  0*5  >  a  good  approximation  to  F  is  F  =  1  +  mv/2,  as 
is  readily  verifieu  from  Eq.  (2),  and  is  also  easily  deduced  from  a  quasi¬ 
static  argument. 

The  asymptotic  limit  of  three  is  quite  unexpected.  To  gain  insight  into 
the  physical  processes,  we  show  in  Fig.  2  the  evolution  of  F  as  a  function  of 
the  normalized  time  t/Tt,  for  various  values  of  v.  It  is  readily  seen  from 
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Fig.  2  that  the  current  enhancement  occurs  mostly  when  Xg  reaches  within  a 
"boundary  layer"  width  of  PR .  The  width  of  this  boundary  layer  (skin  depth) 
is  on  the  order  of  L/v,  as  suggested  by  Fig.  2.  This  may  also  be  deduced  by  a 
dimensional  argument  from  the  diffusion  equation  (1). 

The  above  phenomenon  of  current  enhancement  is  also  observed  when  we 
extend  our  model  to  two  dimensions  (e.g.,  a  current  carrying  wire,  situated 
inside  a  waveguide  filled  with  conducting  medium,  is  given  a  uniform 
transverse  motion).  Specifically,  the  maximum  enhancement  factor  of  three,  as 
well  as  the  presence  of  a  boundary  layer,  is  again  obtained.  Indeed,  these 
features  are  of  such  a  universal  nature  that  they  are  independent  of  the 
initial  location  of  the  current  wire  (or  current  sheet  for  the  one-dimensional 
problem),  as  long  as  >  >  t^. 

This  asymptotic  limit  of  three  demands  an  explanation.  To  fix  ideas,  we 
return  to  the  one-dimensional  model  described  above.  We  let  the  speed  u  be 
finite  and  the  channel  conductivity  a  be  large  (so  that  v  *  au  is  large).  The 
beam  current  always  contributes  one  unit  to  the  enhancement  factor  F,  so  it 
suffices  to  argue  that  the  induced  current  contributes  twice  the  beam  current 
as  the  current  sheet  approaches  PR. 

One  unit  of  the  induced  current  comes  from  the  persistence  of  the  initial 
magnetic  field  3A(x,0)/3x  which  is  already  present  at  t  -  0.  The  associated 
current  distribution,  which  relaxes  slowly  in  the  limit  Td  >>  is 
represented  by  JA  in  Fig.  3*  lb  is  approximately  gaussian  in  space.  The 
total  current  associated  with  is  approximately  equal  to  Ig  since  little  of 
this  initial  magnetic  energy  is  dissipated  on  the  time  scale  t  if  >  >  t^. 


The  remaining  unit  of  the  induced  current  originates  from  the  motion  of 
current  sheet.  [This  unit  is  independent  of  the  initial  A(x,0),  it  therefore 
exists  even  in  the  hypothetical  case  A(x,0)  -  0.]  For  uniform  motion  of  the 


current  sheet,  Eq.  (1)  admits  a  self-similar  solution  of  the  form  Ac(x,t) 

*  f  [x  -  xg(t)3.  Such  a  solution  Ac  yields  an  induced  current  density  given 
by 


Jc  (x,t) 


nv(x-xn)/L 
b  ; 


X  <  x3 
X  >  xB 


(4) 


which  is  sketched  in  Fig.  3.  In  physical  terms,  the  movement  x0  inductively 
drives  a  return  current  at  the  leading  edge  of  Xg,  as  given  by  Eq.  (4). 
However,  Ac(x,t),  which  gives  rise  to  JQ,  does  not  satisfy  the  boundary 
condition  A-(L,t)  -  0.  An  image  current  source  of  the  form  -I.  <5  lx-[2L-x_(t)  ] } 
must  be  added  to  render  the  vector  potential  zero  at  x  -  L.^  Associated  with 
this  image  is  the  current  density 


JD(x,t) 


(1^1)  e  ™(x  -  2L  ♦  x0)/L> 


x  >  2L  -  xB 
x  <  2L  -  xB 


(5) 


which  is  also  sketched  in  Fig.  3.  Summing  up  all  contributions  ( JA, Jg, Jq , JD) 
between  P^  and  PR,  a  current  enhancement  factor  of  3  i3  then  obtained  as 

Xg  ->  L,  since  the  contribution  from  Jq  becomes  negligibly  small,  while  the 

contribution  from  JD  becomes  significant.  From  Eqs.  (4)  and  (5),  as  well  as 
Fig.  3.  it  is  clear  why  the  current  enhancement  would  occur  within  a  boundary 
layer  of  width  L/v  of  PR  when  v  is  large,  as  shown  in  Fig.  2. 

The  above  result  has  been  based  on  idealized  models.  Never theless ,  the 

exposition  of  the  basic  physical  processes  enables  us  to  give  an  assessment  of 
more  complicated  situations.  For  example,  if  the  electrical  conductivity 
a  is  highly  concentrated  only  near  the  central  region  of  the  channel,  (i.e., 


a  is  small  near  the  boundaries  PH  and  P.  ) ,  then  only  Jfl  and  JR  would  be  the 


%  V'  A' ” 


dominant  contributions  as  the  current  sheet  (source)  approaches  P^.  In  this 

case,  an  upper  bound  of  F  -*•  2  is  expected.  As  another  example,  if  the  motion 

of  the  current  source  is  nonuniform,  e.g.,  the  current  source  is  accelerated 

(decelerated)  toward  the  wall  PR,  a  simple  argument  shows  that  the  limiting 

value  of  F  (Fig.  1)  will  be  greater  (smaller)  than  3.  We  have  also  repeated 

-12  2 

our  calculations  with  the  displacement  current  term  (Uirc)  9  A/9t  included  on 

the  left-hand  side  of  Eq.  (1).  We  found  that  inclusion  of  the  displacement 

current  increases  the  enhancement  factor  F  in  general,  but  F  does  not  exceed  a 

value  of  3-2  if  u/c  £  0.3  for  the  simplified  model  treated  above. 

Finally,  we  remark  that  Fig.  2  suggests  the  possibility  of  building  a 

short  pulse  current  amplifier,  the  maximum  amplification  being  a  factor  of  the 

order  of  three.  From  this  figure,  the  pulse  duration  Tp  is  on  the  order  of 

2  2 

t  »  (1/v)  x  L/2u  »  c  /8u  a.  The  transverse  motion  of  the  current  may  be 
P 

provided  mechanically  or  explosively.  This  process  may  also  be  viewed  as  one 
form  of  magnetic  compression. 

In  conclusion,  current  enhancement  takes  place  naturally  when  a  current- 
carrying  body  undergoes  transverse  motion  in  a  conducting  medium.  Especially 
large  current  enhancement  may  occur  when  this  current  moves  to  within  a 
boundary  layer  thickness  from  the  container.  The  recent  experimental 
observations  cited  above  may  be  just  one  of  the  many  manifestations  of  this 
electromagnetic  phenomenon.  Indeed,  because  of  the  simplicity  but  surprising 
nature  of  these  phenomena,  one  may  well  wonder  whether  this  and  related 
phenomena,  in  particular  the  limiting  value  of  F  =  3,  have  been  noted 
elsewhere  (e.g.,  geomagnetic  dynamo,  magnetohydrodynamic  boundary  layers, 
etc. ) . 
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touch  the  conducting  boundary. 
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